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In models where the suppression of the tree level couplings of fermions with Goldstone bosons 
is not protected by symmetry, the loop-induced couplings can in general dominate over the tree 
couplings. We demonstrate this by calculating the decay v — > v' + Majoron in the simplest 
singlet Majoron model. It is shown that 1-loop contributions can be orders of magnitude bigger 
than the tree-level result, and hence allows a wide range of interesting values of neutrino mass 
and mixing to be compatible with cosmological constraints. 



Theoretical models with B — L as a spontaneously broken global symmetry have a Nambu-Goldstone 
boson called Majoron (J), where B and L are respectively the baryon and lepton number. In these models 
a neutrino can decay to a Majoron, v — > v' + J. Such an invisible decay can be of cosmological interest. 
In particular, a neutrino with a mass lying in the cosmologically excluded region between 30 eV and 2 
GeV for stable neutrinos Q would still be allowed if it decays sufficiently fast. In that case the standard 
cosmological constraint derived from the mass density of the universe is B 



m u ^T„/t <30eV, (1) 

where m u and r v are the mass and the lifetime of u, whereas to ~ 5 x 10 17 s is the present age of the 
universe. The bound can be strengthened to M 



m v ^T u /l&s < lMeV, (2) 

provided one assumes the canonical picture of galaxy formation. However, since the physics of galaxy 
formation is not fully understood, this bound should be taken less seriously than the previous one. 

In many simple versions of Majoron Models the decay process occurs with a negligibly small rate ||| 
at the tree level. This follows because to a first approximation the neutrino-Majoron coupling matrix 
often turns out to be proportional to the neutrino mass matrix ||. Thus, in the mass eigenstate basis the 
neutrino-Majoron interaction is diagonal and the decay wouldn't occur at this level. To have a nonzero 
tree- level off-diagonal coupling it is necessary to consider non-leading terms, which are related to the 
mixings of the light neutrinos v and v 1 with some other very heavy species and hence highly suppressed. 

This feature is generally unprotected by symmetries, however. Modifications introduced by loop cor- 
rections can therefore be potentially significant, leading to a rare situation in which the one loop correction 
is larger than the tree-level result. Consider the earliest Majoron model of Chikashige, Mohapatra and 
Peccei § as an example. The model is simple in the sense that it requires only minor modifications over 
the standard model. Therefore the calculation is relatively simple, and the result is dramatic. The model 
differs from the standard model by a right-handed neutrino Nr for each generation, and a Higgs singlet 

$(x) = -^(v 1 + s(x)+iJ(x)), (3) 

which carries 2 units of B — L quantum number. The vacuum expectation value of v%, breaks B — L 
spontaneously and gives rise to a Majoron J. As a gauge singlet, J does not interact directly with the Z 
boson, and hence unlike many other later versions this model has not been ruled out by the measurements 
of Z- width. 

At the tree level, the majoronic decay of neutrino is determined by the Yukawa interaction 

fa _ i 

- £y = V —V aA iP La <j)N RA + V -—MabN la ®N rb + h.c, (4) 

where ipi is the usual lepton doublet and N = C^qN* . <j> is the standard Higgs doublet with its neutral 
component (jp given by 

4>\ X ) = ±=( V2 + H(x)+Mx)). (5) 
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T> and Ai are coupling matrices in flavor space, vi is the vacuum expectation value which breaks the gauge 
symmetry. The mass terms of the neutrinos are given by 



2 ( V 1 Xl 



V 
V T M 



N R 



+ h.c, 



(6) 



where each entry in the mass matrix represents a 3 x 3 matrix for 3 generations of fermions. 

In the present model the elements of T> arise from the same physics that gives rise to the masses of 
charged leptons. As a result, we assume 

VaA < gv 2 Va, A. (7) 

where g is the SU(2) gauge coupling. To proceed, we further assume that CP is conserved and that 

r = v 2 /v! < 1 . (8) 

This implements the well-known see-saw mechanism j?J for the generation of small neutrino masses. In 
addition, we assume for simplicity that 



Mab ~ vi \/A,B. 



(9) 



These assumptions are not crucial. However they do simplify the calculation. 

Calling the doublet Yukawa couplings collectively as y, we can then diagonalize the mass matrix per- 
turbatively in powers of r and y. To order r 2 y 2 one obtains the following relations between the mass 
eigenstates and the weak eigenstates || : 
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ry. To order p , the mass terms are 
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with diagonal mass matrices 
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(10) 
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pMp T ) f, 
F' ( M + \p T pM + ^Mp T p ) F. 



(11) 

(12) 
(13) 



The diagonalization into a light block and a heavy block is performed by the first matrix on the right hand 
side of Eq. (|io|) . The two unitary matrices / and F determine the mixings among the light-light and the 
heavy-heavy neutrinos, respectively. Notice / and F can be complex, even if CP is conserved, unless the 
Majorana eigenstates all have the same CP eigenvalue ||. 

From Eqs. (|j), ( JlOj ) and (12) one can easily show that, to 0(p 2 ) ~ 0(r 2 y 2 ), the xxJ coupling matrix 
is proportional to m which is diagonal ||. Although the next order terms are not diagonal [ 10 1 , they are 
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0(p 4 ) ~ 0(r 4 y 4 ), rendering their contributions to the neutrino decay too small to be of practical interest. 
We shall show that this feature is modified drastically by radiative corrections. To discuss this effect, let 
us parameterize the effective interaction for the decay Xa — ^ Xb + J, keeping terms upto 0(r 2 y 2 ), as 

Asff = -VCb {j^lsha + % a R - T h * a L^j XaJ + h.c, (14) 

in which the first term is the tree-level result and T ba comes from loop contributions to be calculated below. 
The 1-loop diagrams for the decay in the Feynman gauge are shown in Figs, [j], and [3|. Graphs with an 
H-s mixing are negligible due to (Mj. Charged-current interactions do not make any contribution at this 
level. The Feynman rules for relevant couplings are summarized in Fig. [|. 

The contributions of various diagrams are organized in accordance with the virtual bosons propagating 
in the loop. For example, T^ z > represents the contribution from all diagrams containing a virtual Z. 
Neglecting the mass of the final state neutrino and using the shorthanded 

P = f T VM~ l F* , (15) 

we find, upto terms of order 0(r 2 y 2 ) times a potential logarithmic factor, 

T {Z,H, V ) = 0; (16) 

Tt ] = ^E^^(^) 3 /(Mi,M s 2 ), (17) 
^ = ^En^(^) 3 /(Mi,0), (18) 

T ba ~ ^l^P^PaA^ — J M 2_ M 2^M2> ( 19 ) 

where M s is the mass of the singlet scalar s and 

f(j '» )= f * «+-)«+») (20) 



which is divergent. However, the sum of Eqs. ( p^6[]19| ) is finite: 

T ba = 2T b [ Js \ (21) 

Using the definition of P from Eq. ([lq), one sees that Tba is of order m a /Ait 2 v\ times an interfamily 
mixing for Ma ~ V\ and reasonable choices of M s /Ma- Despite the loop suppression factor 1/4-7T 2 , for a 
wide range of interesting values of m and M this result is orders of magnitude bigger than the tree-level 
result, which is C(r 4 y 4 ), as pointed out earlier. 

There is a technical subtlety associated with the choice of the representation of Majoron fields. Instead 
of using Eq. (H|) as the definition of J, one can use a non-linear representation: 

*(x) = 4=(«i + s{x))e iJ ^' Vl . (22) 
v2 

The interaction Hamiltonian now has momentum-dependent terms. As a result, the formulation of Feyn- 
man rules is different Jll]] from the usual situation. The two representations, Eqs. @ and (^2|), give 
identical results for tree-level processes in which momentum-dependent interactions are negligible. For the 
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case at hand, a naive calculation following ( p2|) without taking into account the associated Feynman rule 

f Js) 

difference would incorrectly give T ba = 0, and as a result the final answer would be a factor of 2 smaller. 

Let us now examine the significance of the loop result. With our assumption of CP invariance, the 
mixing matrices satisfy the constraints || |l2| f* b = r?b/ a 6, and a similar one for Fab, where the CP 
eigenvalue of Xb is with r\ b = ±1. Thus, 

Pl A = r ]a r ]A P aA . (23) 

Due to the Majorana nature of the neutrinos, the hermitian conjugate term in Eq. (|14|) contributes equally 
as the other one. The amplitude of the decay can therefore be written as 

*4(Xa ->• Xb + J) = ^0 ba u b (R - r] a r] h L)u a , (24) 
where u aj b are the standard spinors, and for convenience we have defined a parameter O^a'. 

which plays the role of m a jv\ times an effective light-light interfamily mixing. Note that 9 vanishes when 



vi — > oo, since the heavy neutrinos decouple in that limit. Eq. ( |24| ) demonstrates that the transition 
amplitude is purely pseudoscalar if the two Majorana states have the same CP property, and purely scalar 
if their CP eigenvalues are opposite. Intuitively, this is obvious since the CP eigenvalue of the Majoron is 
negative. 



From Eq. (24), the rate of the decay for either sign of ?y a r/b is 

Ohm. 



Substituting ( ^6|) into (Q) yields 



r-\xa^XbJ) = ^- (26) 



9 2 ba (leV/m a ) > 1.5 x 1(T 32 . (27) 



This result implies that, for a wide range of v\ and neutrino masses and mixings, the radiatively-induced 
neutrino-Majoron decay is fast enough to allow the neutrino be compatible with the big bang cosmology. 

To have a feeling for our result, we take from (0), (|TJ) and (||) our order-of-magnitude estimates 
m a ~ v\p 2 and 9 ~ p 2 . Replacing p by rap/vi, where TO£> is the relevant Dirac mass, we obtain m a ~ 
10 4 eV x (1 TeV/^i) x (m^/100 MeV) 2 . Thus for v\ = 1 TeV and = m^, m a ~ lOkeV. Such a neutrino 
can easily satisfy the cosmological constraint (||) which only requires from (27) 

m a > 1.5- 10" 8 eV x (^i/lTeV) 2 . (28) 

For vi of order TeV, even after taking all the uncertainties of an order-of-magnitude estimate, this shows 
that the mechanism can provide a fast decay for neutrinos of any mass. One can also express Eq. (28) as 
a bound on vi, i.e., v\ < 10 4 TeV x (m/j/100 MeV) 2//3 . If one uses the bound of galaxy formation from Eq. 
(|2[), the constraint is similar to Eq. ( |27| ) with the right hand side replaced by 6.6 x 10 -28 . The corresponding 
bound on the neutrino mass is m a > 6.6-10 -4 eV x (vi/1 TeV) 2 with the same phenomenological conclusion. 

On the other hand, if one has to rely on the tree level contribution only [Jl(J, instead of 8/2tt 2 we 
obtain 2p 4 as the order-of-magnitude estimate for the coefficient in the amplitude of Eq. (|24|), the factor 
2 occuring due to the Majorana nature of the neutrinos. The bound in Eq. (|l]) can then be satisfied if 

m a > 20keV x ( Vl /1 TeV) 4/3 . (29) 



This can be satisfied for a narrow range of values of the mass of z/„ provided v\ is not much larger than a 
TeV. However, for the bound in Eq. (0), the factor 20 keV from this equation should be replaced by 750 
keV, which cannot be satisfied for any j/„ mass allowed by experiments. 

In conclusion, we have calculated 1-loop radiative corrections to the neutrino-Majoron decay in a 
singlet Majoron model. The result shows that loop contributions can decisively dominate the tree-level 
contribution, allowing a wide range of interesting values of neutrino mass and mixing be compatible with 
the big bang cosmology. Note that the sole reason that the 1-loop result dominates over the tree level 
result is just because the tree level suppression is accidental in nature, that is, it is not enforced by any 
symmetry. It has little to do with the other details of the model. Therefore, one expects that the same 
conclusion can be drawn in a large class of models in which there is such accidental tree level suppression 
of off-diagonal majoron couplings. In addition, since the accidental suppression disappears at the 1-loop, 
one does not expect that including higher loop contributions would change our result qualitatively. 
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Figure Captions 



Figure 1: Onedoop graphs for the leading terms of the decay Xa —> Xb + J involving a Z, (p or H internal 
line. The contribution of these diagrams appear in Eq. (|l6|). 

Figure 2: Same as Fig. |l| except the internal boson lines are either s or J. The contribution of these 
diagrams are summarized in Eqs. ( |l~7l) and (|i~8|). 



Figure 3: Diagrams involving the trilinear scalar vertex JJs which contribute at the leading order. Their 



contribution is summarized in Eq. (19). 



Figure 4: Feynman rules relevant for the calculation of ldoop corrections to the decay Xa ^ Xb + J ■ The 
Feynman rules for fermion vertices with H (s) can be obtained by multiplying the corresponding <p (J) 
vertex by —475. 
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